Inhomogeneous magnetism in the doped kagome lattice of LaCuO 
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The hole-doped kagome lattice of Cu 2+ ions in LaCu02.66 was investigated by nuclear quadrupole 
resonance (NQR), electron spin resonance (ESR), electrical resistivity, bulk magnetization and spe- 
cific heat measurements. For temperatures above ~ 180 K, the spin and charge properties show an 
activated behavior suggestive of a narrow-gap semiconductor. At lower temperatures, the results 
indicate an insulating ground state which may or may not be charge ordered. While the frustrated 
spins in remaining patches of the original kagome lattice might not be directly detected here, the 
observation of coexisting non-magnetic sites, free spins and frozen moments reveals an intrinsically 
inhomogeneous magnetism. Numerical simulations of a 1/3-diluted kagome lattice rationalize this 
magnetic state in terms of a heterogeneous distribution of cluster sizes and morphologies near the 
site-percolation threshold. 

PACS numbers: 



I. INTRODUCTION 

Whether the physics of quantum spin liquids is relevant 
to the occurrence of high temperature superconductivity 
in copper oxides has been one of the most debated issues 
in condensed matter physics over the last two decades. 
This question has triggered important research efforts 
in two directions: i) the physical realizations of quan- 
tum spin liquids, and ii) the doping of such systems with 
charge carriers. Natural candidates for realizing super- 
conductivity from a spin liquid are compounds such as 
transition-metal oxides where unpaired electrons reside 
on dimer, ladder or triangle-based networks. 1 In general, 
two-dimensional lattices are more favorable than ID sys- 
tems for which localization is expected at low tempera- 
ture (T). 




FIG. 1: Copper (large light spheres) and oxygen (small dark 
spheres) sites in the CuOo.66 kagome planes of LaCu02.66-~ 
These planes alternate with layers of LaOe octahedra. 



Copper oxides RCu0 2+(5 (R = Y, La, In, Sc, Nd^ 
with a delafossite structure^ provide an interesting play- 
ground for studying the effect of doping in a 2D lattice 
where Cu ions form triangles. Depending on the oxygen 
content S, a variety of triangle-based magnetic networks 
can be obtained as the additional oxygen O^ occupies 
different positions in the Cu planes (i.e. the geometry 
of the magnetic network and the formal Cu valence are 
necessarily linked together). Substantial progress in the 
synthesis of these materials has recently led to the iden- 
tification of well-defined structures^ 

Among copper dclafossites, LaCu02.66 appears to be 
of particular interest. It consists of layers of edge-sharing 
LaOg octahedra alternating with triangular planes of 
Cu 2+ ions. With respect to LaCu02, the extra oxy- 
gen ions are located in the Cu planes. For the con- 
centration S = 2/3, the oxygen ordering is such that 
the magnetic network connecting the Cu ions has a per- 
fect kagome geometry, made of corner sharing equilat- 
eral triangles with an oxygen at their center (Fig. [!])£ 
This ideal situation should be contrasted with that in 
the almost isostructural compound YC11O2.66 which will 
be mentioned several times in this paper. The CuOo.66 
planes of this yttrium counterpart also form a kagome 
lattice but the distorted structure probably leads to non- 
negligible next- nearest neighbor (i.e. inter-triangle) ex- 
change interactions^ The highly frustrated kagome lat- 
tice, of which there exists a number of good, but never 
perfect, experimental realizations (see ie{s£~— for just a 
few examples) is one of the most popular candidate for a 
quantum spin liquid in two dimensions. 

With a formal valence of Cu ions of +2.33, LaCu02.66 
(as well as YCuC>2.66) should be considered either as a 
kagome lattice doped with 33 % of itinerant holes or 
as a diluted kagome lattice with 33 % of non-magnetic 



(Cu 3+ ) sites (Note that band structure calculations indi- 
cate that the states at the Fermi level have predominantly 
Cu d z 2_ r 2 character—). Both possibilities are of consid- 
erable theoretical interestj 11 ' 12 Although LaCu02.66 has 
received early theoretical attention ; 10 ^ 3 little is known 
concerning its electronic properties. For 8 = 0.64, which 
can be viewed as a disordered version (i.e. with 0.02 
oxygen vacancies) of the stoichiometric 5 = 0.66, early 
studies found a metallic ground state with strong cor- 
relations and strong disorder effects close to the metal- 
insulator transition s 14 ^ 5 For the more recently synthe- 
sized d = 0.66 concentration, which presents improved 
structural homogeneity, a muon spin rotation (/iSR) 
study found antiferromagnetic (AFM) order in ^20% of 
the sample volume below 40-50 K.— 

The paper is organized as follows: Sections II, III and 
IV arc devoted to experimental results in LaCu02.66- 
Section V presents numerical simulations of the kagome 
lattice at a 1/3 dilution level, a description justified 
by the experimental evidence that the doped holes are 
strongly localized at low T. Section VI summarizes the 
results. Key issues and open questions are discussed in 
section VII. 



II. EXPERIMENTAL DETAILS 

Polycrystalline samples were used for all measure- 
ments. They were prepared by solid state reaction fol- 
lowed by an oxidization procedure with long annealing 
times at the lowest possible temperatures in order to get 
good quality samples.— The results from magnetization 
and specific heat measurements showed minor sample de- 
pendence at the quantitative level, presumably due to the 
different amount of impurities, but not at the qualitative 
level. The magnetization data shown in this article were 
obtained on the sample used for NQR measurements and 
for a ^SR studyJ£ 

The specific heat was measured by the pulsed relax- 
ation method and the resistivity was measured by the 
standard four-contact method, both in a commercial 
setup (Quantum Design PPMS). The magnetization was 
measured by the axial extraction method in a home-built 
magnetometer in the temperature range 2-800 K and in 
magnetic fields up to 10 T as well as on a Quantum De- 
sign SQUID magnetometer (MPMS) in the range 2-300 K 
and in fields up to 5 T. 

The ESR measurements were performed at different 
frequencies and temperatures using a home built multi- 
frequency high-field spectrometer. Different frequencies 
were investigated (54—115 GHz) using back- wave oscilla- 
tors and Gunn diode sources, guided to the sample with 
oversized waveguides in a transmission probe. The spec- 
tra were recorded as a function of magnetic field at a 
fixed frequency. 

The 139 La and 63,65 Cu NQR measurements were per- 
formed on a home-built pulsed spectrometer. Let us re- 
call that NQR measurements, unlike nuclear magnetic 
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FIG. 2: (Color online) Main panel: electrical resistivity vs. 
temperature. The line represents a thermally activated be- 
havior (see text). (Inset) Same data vs. inverse temperature. 



resonance (NMR), are performed in zero external mag- 
netic field. The signal intensity was determined from 
the integral over time of the spin-echo signal, and di- 
vided by f 2 as a function of NQR frequency /. The 
spin-lattice relaxation time (Jf) was estimated from a 
single-exponential fit of the time dependence of the nu- 
clear magnetization M z after a single saturation pulse: 
M z {t) = M(t = oo)(l - exp(-(3£/Ti) a ))J£ The stretch- 
ing exponent a was found to vary smoothly from 0.9 at 
300 K to 0.7 at 4.2 K, indicating a moderate distribution 
of Ti values.i&^o 



III. MACROSCOPIC MEASUREMENTS 



A. Resistivity 

The resistance rapidly increases on cooling and be- 
comes impossible to measure reliably below ^150 K, most 
probably because of a large contact resistance. The in- 
trinsic transport properties of this oxide are thus diffi- 
cult to determine from bulk measurements on a pressed 
pellet. As previously reported in LaCu02.64^, the resis- 
tance of LaCu02.66 does not follow a simple activated 
behavior over the whole temperature range (Fig. [2]). A 
fit to an activated law p — po cxp(An/ ksT) can nonethe- 
less be performed over limited temperature intervals be- 
tween 300 and 180 K. This yields an activation energy 
A p = 2700 ± 300 K. Although the T dependence may be 
affected or even dominated by the grain boundaries, an 
activated law is also found in other measurements (see 
below) and is thus consistent with a semi-conducting be- 
havior. 



B. Bulk magnetization 

The bulk magnetization M was found temperature re- 
versible up to 750 K, above which the sample is subject 
to chemical transformation. M/H is equal to the linear 
susceptibility \ down to 5 K, below which the isother- 
mal magnetization versus field starts to deviate from lin- 
earity. As can be seen in Fig. [3J M/H displays a non 
monotonic temperature dependence. The T x M/H ver- 
sus T representation of Fig. [4] gives another view on the 
different regimes. The data suggest the presence of two 
contributions with opposite T evolutions. Starting from 
high temperatures, the decrease of M/H on cooling is 
tentatively modeled by an activated behavior. Then, the 
sharp rise of M/H below ^200 K can be accounted for by 
a Curie- Weiss law. The whole susceptibility then writes: 



X = Xo + C/(T + 9) + a x exp(- 



A x /k B T), 



(1) 



and the following values are obtained from fitting of 
the data to this expression: a constant term due to 
sample environment and diamagnetic contribution xo = 
—4.5 x 1CP 6 cmu/(mol.Oc), the Curie constant C = 
0.016 ± 0.002 emu.K/(mol.Oe), the Curie- Weiss temper- 
ature 9 = 13.8±5K, o= 1.34 x 10" 3 cmu/(mol.Oc), and 
the susceptibility gap A x /k B = 1281 ± 200 K. 

If a gap were present in the spin excitations, its max- 
imum possible magnitude would be set by the main 
magnetic exchange coupling J. The intra-triangle J in 
LaCu02.66 is not known but it is expected to be of the 
order of a few hundred Kelvin, since J values range from 
170 to 660 K in YCUO2.5, a complex two dimensional 
array of similar Cu 2+ triangles with an oxygen at their 
center i 21 i 22 Inter-plane magnetic couplings, on the other 
hand, are considered to be much weaker than in-plane 
couplings, despite the d z 2_ r 2 orbital occupation.— It is 
thus unlikely that the large gap magnitude A x ~ 10 3 K 
can be explained by a spin gap. Instead, the similarity in 
the orders of magnitude of the gap values suggests that 
the activated behaviors in the resistivity and suscepti- 
bility have a common origin. The simplest scenario we 
can think of is a thermal activation of states in the va- 
lence band of a small-gap semiconductor (see for instance 
ref. 23 ). The rise of % at high temperatures then results 
from the increasing contribution of the Pauli susceptibil- 
ity of the conduction electrons. 

The magnitude of the Curie constant C indicates that 
the Curie- Weiss behavior of the susceptibility arises from 
5% of spins 1/2 (relative to the total number of Cu ions). 
This term will be shown below to be mostly of intrinsic 
origin. The fitting, however, does not capture fine de- 
tails of the shape of M/H, especially the non linearity 
of the inverse of M/H at low T (inset of Fig. [3J. This 
can be due to the contribution from collective spins 1/2 
from small clusters with a complex geometr y 24 ' 25 and/or 
to the onset of partial magnetic ordering whose existence 
is supported by the slight difference between field-cooling 
and zcro-ficld-cooling below «80 K. This, together with 




FIG. 3: (Color online) M/H vs. temperature T. Data above 
and below 300 K results from two separate measurements: 
in the temperature range [2-300 K], M was measured in a 
field #=10 kOe while in the range [300-700 K], M/H was 
extracted from linear fits of M versus H up to 60 kOe. The 
two methods give the same value of M/H at 300 K. The 
line is a fit following Eqn. (1). Dashes represent the function 
a x exp(— A x /fc_eT) with A=1281 K. The inverse curves are 
displayed in the inset. 
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FIG. 4: (Color online) (a): M/H times T versus T, in 
H=W kOe. (b) M/H vs. T measured in 100 Oe after zero- 
field cooling (open circles) and after cooling in a field of 100 
Oe (filled circles). Open and filled squares are the same data 
from which the Curie- Weiss contribution xo + C /(T + 6) has 
been subtracted (see text). The residual upturn below fa 13 K 
is possibly due to paramagnetic impurities. 



a progressive downturn of T x M/H on cooling (Fig. @| 
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FIG. 5: (Color online) specific heat C p on a log-log scale at 
different magnetic field values. Solid lines are fits of the data 
(symbols) withC p = fiT 3 + /3 2 T 2 + aR(-^) 2 exp-(^)/(l + 

e~~&f with /3=0.00015 J.K.moP 1 , /3 2 =0.00294 J.moP 1 , 
«=0.0064 J/mole K 2 and the values of A are shown in the 
inset as a function of magnetic field. The dashed line shows 
the T 3 + T 2 contribution. 



of various spin liquid phases of the 5 = 1/2 Heiscnbcrg 
model on the kagome lattice. 28 ! 29 It could also arise from 
excitations of the 2D-AFM order. 

- a Schottky anomaly, likely due to the presence of dis- 
crete low energy levels, accounts for the field dependence 
of C p . The deduced energy barrier varies linearly with 
the field (inset to Fig. [5]). Assuming that this terms arises 
only from two-level paramagnetic impurities, it would ac- 
count for 0.6 % of the magnetic Cu. However, the ob- 
servation of a residual contribution in zero field reveals 
that the Schottky term docs not solely originate from im- 
purities but also from the distribution of discrete energy 
levels characterizing the small clusters. 

- a temperature-linear term was tentatively introduced. 
However, we could clearly better fit our data without this 
last term. Since this term is attributed to conduction 
electrons in LaCu02.64^, its absence here is another in- 
dication that LaCu02.66 nas an insulating ground-state. 



IV. LOCAL PROBE MEASUREMENTS 



'La NQR 



indicates that the viscosity of the spin correlations pro- 
gressively increases due to some sort of freezing in the 
spin system. However, the relatively modest amplitude 
of the effects observed in the susceptibility (in particular 
as compared to YCu02.6a^) is inconsistent with a con- 
ventional transition towards bulk long range order. This 
observation agrees with a spin freezing in only ~20% of 
the sample volume as inferred from //SRJ£ 



C. Specific heat 

At high temperature, the specific heat C p measured in 
zero magnetic field approaches the expected Dulong and 
Petit limit. No sharp transition is detected in the investi- 
gated temperature range, thereby ruling out the presence 
of a conventional, homogeneous, long-range magnetic or- 
der (sec Fig. [5]). A field dependence is observed below 
10 K. No evidence for an activated behavior, that could 
reveal a large spin gap, was found in these data. The 
data below 11 K were fitted by the sum of several con- 
tributions (see Fig. [5}: 

- a T 3 term accounting for the Debye lattice, which 
cannot be evaluated independently here because there is 
no non-magnetic isostructural reference compound. 

- a T 2 term, as previously reported for LaCu02.64^ 
The T 2 term in LaCu02.66 is most probably of mag- 
netic origin since we found it to be absent in the specific 
heat of LaCuC>2 which, although not isostructural with 
LaCuC>2.66j presents a similar lamellar structure. Such 
a T 2 law has been observed for instance in the kagome 
bilayer of SCGO even at large dilution levels 2 ^ and has re- 
ceived theoretical interpretations in terms of excitations 



Starting from its low frequency side, the NQR spec- 
trum of LaCu0 2 .66 at T = 4.2 K (Fig. |5J) first shows 
that there are three different 139 La sites from the elec- 
tric field gradient point of view. The relative intensities 
of these Lai, La2 and La3 sites represent approximately 
27, 13 and 60 % of the total La intensity. In princi- 
ple, there arc six different crystallographic sites for La 3+ 
ions in this material^ Nevertheless, these six sites can be 
grouped into three categories, according to the number 
(six, seven or eight) of oxygen anions to which each La 3+ 
is coordinated^ The relative population of these three 
types of coordination cannot be accurately determined 
because of uncertainties regarding the complex stacking 
sequence of copper-oxygen planes and the occurrence of 
numerous stacking faults.— Still, assuming an equal prob- 
ability for the three possible stacking sequences (named 
A, B and C in Ref£), the sixfold, sevenfold and eight- 
fold coordinated La would represent 12, 44 and 44 % 
of the sites, respectively. Given the large uncertainties, 
these numbers appear compatible with the relative NQR 
intensities so that the three different 139 La NQR sites 
could indeed correspond to the three different coordi- 
nation numbers of La 3+ ions. It cannot be excluded, 
however, that part of the differentiation of electric field 
gradients at the La site is related to the distribution of 
charges in CuO^ planes. 



B. 63 ' 65 Cul NQR 

Around 16-17 MHz, a single, intense, NQR line is ob- 
served for each of the two isotopes 63 Cu and 65 Cu. This 
NQR site is named Cul in Fig. [5] From their 63 Cu NMR 
line in LaCu02.64, Walstedt et al. deduced a quadrupolc 
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FIG. 6: (Color online) NQR spectrum of LaCuC>2.66 at T — 4.2 K. (a) Three different 139 La sites (grey peaks) and a main 
63,65 Cul site are resolved. The three La sites are observed in two different sets of three lines, that correspond to the ±| -H- ±| 



and ±| <-> ±| 



transitions within the energy levels of the nuclear spin 139 J = \. 



The ±| •(-> ±| transitions occur at too low 



frequency to be observed with our setup. Both Cu and Cu have a spin / = 3/2, thus only one ±-| o § NQR transition for 
each isotope/site, (b) Cu2a and Cu2b sites corresponding to minority paramagnetic sites (see text). The total intensity of these 
sites represents a few % of the Cul intensity, (c) At T — 1.6 K, the broad ' 5 Cu resonance around 121 MHz is attributed to 
sites in presence of a 10.2 T internal field, i.e. magnetic order. The absence of data points around the peak maximum is due 
to the presence of a spurious signal from the NQR probe. 



coupling constant vq = 15.5 MHz and an asymmetry pa- 
rameter of the electric field gradient tensor r\ = 0.65 at 
T = 293 K£ These parameters imply a NQR frequency 
^nqr = vq^jV + 77 2 /3 = 17.7 MHz quite close to our 
measured value of 17.55 MHz at 300 K (Fig. W)- This 
suggests that the Cul site in LaCu02.66 could have a 
similar electronic environment at high temperature as the 
site observed in LaCu02.64- 

The frequency of the Cul line is constant a low T but 
it acquires a marked temperature dependence (best de- 
scribed as linear) above 150 K (Fig. [7K). The date are 
consistent with an electronic crossover near 150 K but a 
continuous variation throughout the temperature range 
(such as an activated T dependence shown in Fig. [3t 
with a characteristic energy E ~ 620 K) cannot be fully 
excluded. 

The spin-lattice relaxation rate 1/Ti of LaCu02.66 has 
similar values at 300 K as in LaCu02.64(Fig. [7]). On 
cooling, however, the two compounds show a very differ- 
ent temperature dependence. In LaCu02.64, 1/Ti values 
remain high down to low Tji£ consistent with a metal- 
lic ground state. In LaCu02.66j on the other hand, the 
particularly long T\ of 14 s at T = 4.2 K reveals that 
the Cul site is non-magnetic. This site should thus cor- 
respond to the non-magnetic Cu 3+ ions expected from 
the localization of holes at low temperature. In princi- 
ple, Cu 2+ spins in dimer arrangements could also have a 
very small 1/Ti at low T but they should have a different 
resonance frequency. 

Above 150 K, the strong temperature dependence of 
1/Ti can be fit to 1/Ti oc cxp(-A T i/fc B T) with A T1 ~ 
1216±14K (Fig.[7b,c). Within error bars, this gap value 
is equal to that inferred from the magnetic susceptibility, 
while being somewhat larger than that inferred from the 



resistivity (let us recall that the value of this latter might 
be affected by the grain boundaries). Since there is no 
doubt that the NMR signal comes from the bulk of the 
sample and is thus intrinsic, the comparable gap magni- 
tudes tend to suggest that the activated dependence of 
the resistivity measured in this temperature range could 
indeed be intrinsic and not solely due to grain bound- 
aries. Below 150 K, the deviation from the activated be- 
havior indicates that the relaxation becomes dominated 
by residual relaxation processes (such as relaxation by 
paramagnetic impurities, quadrupole relaxation, etc.). 



C. 63 65 Cu2a,b NQR 

In the 20-30 MHz range, two other Cu NQR sites 
(Cu2a and Cu2b) are found. Their intensity integrated 
over frequencies is weak, at most a few pcrcents of the 
Cul intensity. At 4.2 K, the spin-lattice relaxation time 
T\ of Cu2b is 24 ms, that is shorter than that of Cul by 
almost three orders of magnitude, indicating that, unlike 
Cul, these sites carry a finite spin polarization. In agree- 
ment with the above analysis of the Curie- Weiss contri- 
bution to the spin susceptibility, these Cu2a and Cu2b 
sites should correspond to intrinsic paramagnetic spins 
1/2 in the sample, rather than to diluted defects. 



D. 



5 Cu zero-field NMR signal 



At T = 1.6 K and at frequencies around /o ~121 MHz, 
a broad signal is attributed to Cu nuclei in the presence 
of an internal, static, magnetic field H mt = 27r/o/ 63 7 = 
10.2 T. Such a value of the resonance frequency is typical 
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FIG. 7: (Color online) (a) T dependence of the NQR fre- 
quency of the Cul site, with a clear saturation below 150 K. 
The continuous (blue) line is a linear fit to the data from 
150 K to 300 K. The dashed (black) line is an exponential fit 
a+bexp(-E/k B T) with E = 622 K. (b) T dependence of l/7\ 
for the Cul site, compared to data in LaCu02.64 from ref.— . 
The solid line is a fit to an activated behavior exp(— Aii /ksT) 
with A^ = 1216 K, which holds from 300 K down to approx- 
imately 150 K. (c) Same 1/Ti data and fit in log scale vs. 



for AFM order in copper oxides, such as the high tem- 
perature superconductors. This signal is thus a zero-field 
NMR signal rather than a NQR signal since the dominant 
source of splitting of the nuclear levels is the Zceman in- 
teraction rather than the quadrupolar one. The presence 
of frozen moments on some of the Cu sites of the sample 
agrees with a muon spin rotation (/iSR) study performed 
on the same sample,— which indicates that AFM order 
is present below 40-50 K in ^20% of the sample volume. 
Confirmation of this number was not possible from the 
present nuclear resonance data. This signal was found to 
disappear very rapidly on warming (it was much weaker 
at 4.2 K) most probably because of dramatic shortening 
of the transverse relaxation time T^. The NQR signal 
from these sites remains to be observed in the paramag- 
netic phase. 



E. ESR 

The ESR line at 95 GHz is a single Lorentzian, in- 
dicating that any anisotropy of the g factor is not re- 
solved at this frequency. Its intensity at low temperature 
(not shown) follows a Curie- Weiss law suggesting that 
this signal stems from the same spins as probed in the 
bulk susceptibility. As Fig. [5] shows, the temperature 
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FIG. 8: (Color online) Field B a and width AS of the ESR 
resonance line (shown in the inset for T = 5 K) at a frequency 
of 95 GHz. The width is defined as the peak-to-peak distance 
of the signal shown in the inset. Vertical dotted lines show 
temperature crossovers. 



dependence of both the field and the width of the res- 
onance shows a first crossover near 150-180 K followed 
by a sharper change near 30-40 K (for each of these two 
temperature ranges, the lower T value is that indicated 
by the resonance field data, while the larger one is given 
by the width data). The higher temperature crossover 
should be related to the change observed in NQR fre- 
quency in the same temperature range. The additional 
shift /broadening of the resonance field below 30-40 K 
is clearly due to the magnetic freezing found to onset 
around 40-50 K in ^SRpi£ These preliminary ESR re- 
sults demonstrate that the ESR signal is from intrinsic 
paramagnetic spins in LaCuQ2.66- 



V. NUMERICAL SIMULATIONS 

A. Cluster sizes 

The experimental evidence for localization of the 
doped holes at low temperature signifies that LaCu02.66 
should be described as a 1/3-diluted kagome lattice. We 
thus performed a number of numerical simulations in or- 
der to provide a concrete picture of the system at such 
a high dilution level d = 0.33, which turns out to almost 
correspond to the site percolation threshold p c =0.6527 of 
the kagome lattice^ (i.e. d c = l —p c = 0.3473). The dis- 
tribution of magnetic cluster size (Fig. [9]) was obtained 
by averaging over a large number of randomly depleted 
lattices and it was checked to converge rapidly when in- 
creasing the lattice size. At large dilution levels, only 
small clusters (typically less than 10 sites) are present. 
For intermediate dilution d ~ 0.5, there is a homogeneous 
distribution of sizes up to typically 100 sites (see inset of 
Fig. 9). 
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FIG. 9: (Color online) Coverage of the lattice (= Number of 
clusters of a given size times the number of sites in the cluster 
divided by the total number of lattice sites), as a function 
of cluster size (i.e. number of sites in a magnetic cluster), 
for different values of the dilution level d. The inset shows 
an expanded view at low values of cluster sizes. The data 
are from simulations of a planar lattice of 2700 sites. Since 
no interplanar exchange coupling is considered, the lattice 
consists of a single plane. 



As the percolation threshold p c = 1 — d c is crossed, the 
intermediate-size clusters appear to coalesce into large 
clusters and the distribution becomes essentially bimodal 
with a minority of sites in small clusters and a large ma- 
jority of sites embedded in extended clusters, some of 
which are percolating. An example of a diluted kagome 
lattice of N = 1200 sites at d = 0.33 is shown in Fig. [TO] 
At this concentration and for N = 2700, clusters of less 
than 100 sites cover about 12 % of the lattice while clus- 
ters of more than 1000 sites cover 55 %. 



B. Cluster morphologies 

Fig. [TU] shows that, to first approximation, two very 
different types of connectivity can be singled out in the 
magnetic clusters. On the one hand, there are spins for 
which the local environment remains identical to the pris- 
tine kagome lattice. On the other hand, there are spins 
in a weakly frustrated environment, including chain- like 
fragments. We hypothesize here that the weakly frus- 
trated spins will have a tendency towards AFM-like or- 
dering whereas the more frustrated sites will still form 
spin liquid-like arrangements (which may be of various 
kinds: valence bond crystals, dimer liquid, valence bond 
glass, etc.) and thus remain mostly non-magnetic with a 
susceptibility relatively close to that of the pure kagome 
lattice. 

Furthermore, a similar description has already been 
proposed to explain the existence of two different sites in 
the Mossbauer spectra of LuFeMg04, a triangular lattice 
of Fe 3+ spins with a site dilution level of 50 % coinciding 
with the percolation threshold of this latticed The mag- 
netic freezing of this compound (around 33 K, that is well 



below the Curie- Weiss temperature of -500 K) has thus 
been ascribed to the freezing of the weakly frustrated 
spins. Although the ground state of the triangular lattice 
(long-range ordered with spins oriented at 120° from each 
other) is different from that of the kagome lattice (non- 
magnetic with short-ranged spin correlations), there are 
obvious analogies with the situation in LaCu02.66- 

Of course, this description is a caricature of the actual 
situation since there exist many different morphologies 
with all possible intermediate levels of frustration in be- 
tween these two limits. Nonetheless, this simple descrip- 
tion helps to understand qualitatively the inhomogeneous 
magnetism in LaCu02.66 and a more sophisticated model 
is beyond the scope of this paper. Note also, as the next 
section will emphasize, that the local environment and 
the symmetry of sites belonging to large clusters is irrel- 
evant from the point of view of the susceptibility. 



C. Calculation of the magnetic susceptibility 

To first approximation, the Curie-Weiss susceptibil- 
ity arises from isolated 5=1/2 sites and from orphan 
spin o 33 i 34 within clusters with an odd number of sites. 
In fact, it can be shown that the largest clusters (i.e. 
those with the most pronounced kagome character) have 
a negligible contribution to the susceptibility in the tem- 
perature range of interest here as their Curie-like tail 
becomes non-negligible only at very low T. A practical 
cutoff value of 13 was chosen for the size of the clusters 
considered in the calculation of the susceptibility. Each 
of these clusters of N sites (N < 13) was approximated 
to be a chain of N sites, whose magnetic linear suscep- 
tibility was then calculated from exact diagonalization. 
The linear topology for the small clusters is obviously a 
crude approximation. In particular, there is a non linear 
inverse susceptibility yielded by the progressive construc- 
tion of collective spin 1/2 resulting from entanglement in 
small clusters with triangular geometryi 24 i 25 i 36 However 
at low T, the entangled spins 1/2 should yield a Curie- 
Weiss behavior. Finally, an average susceptibility was 
calculated by considering the distribution of small clus- 
ters in a set of randomly depleted kagome lattices, for 
each dilution level. The susceptibility computed with 
these assumptions was fitted with a Curie- Weiss model. 
The resulting percentage of spins 1/2 per total number 
of Cu is shown in Fig. [IT] as a function of d. For 1/3 di- 
lution, the value of ~3 % matches well the experimental 
result of 5 % (which also includes the contribution from 
extrinsic defects) (Section IV. B). 



VI. SUMMARY OF MAIN RESULTS 

A. From experiments 

As far as magnetic properties are concerned, our ex- 
perimental results can be summarized in terms of co- 




FIG. 10: (Color online) Example of cluster distribution from a 33% random site dilution on the kagome lattice with 20 x 20 x 3 = 
1200 planar sites with three sites per unit cell. Examples of a small cluster (red line), and of a non- frustrated environment 
(blue line) and of a frustrated environment (green line) in the same big cluster are shown. Not fully apparent in this figure is 
the fact that clusters of less than 100 sites cover approximately 12 % of the lattice while clusters of more than 1000 sites cover 
~55 %. 
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FIG. 11: (Color online) Calculated percentage of Cu sites 
contributing one spin 1/2 to the Curie- Weiss susceptibility, 
as a function of the dilution level per Cu site d. The vertical 
line at d c = 0.3473 indicates the site percolation threshold 
p c =0.6527 of the kagome lattice (d c = 1 — p c ). See text for 
details. 



existence of different Cu sites at low temperature in 
LaCu0 2 .66: 

(1) Cu 2+ sites experiencing magnetic order. We con- 
jectured that these sites, which are observed in nuclear 
resonance and /xSR experiments^ belong to extended 
parts of the clusters, such as chain fragments, whose ge- 
ometry effectively releases the frustration. 20 % of the 
implanted muons experience a static field from these mo- 
ments, suggesting that the magnetic sites represent a 



fraction of the sites which is of the order of 20 %, ie 

(2) Isolated paramagnetic Cu 2+ spins. These sites can 
account for the Curie-like tail of the susceptibility, possi- 
bly for the Cu2a and Cu2b NQR sites, and for the Schot- 
tky anomaly in the specific heat. We have argued that 
they originate from small clusters, although some of them 
may also be of extrinsic origin. 

(3) Non- magnetic Cu 3+ . All, or part, of the Cul NQR 
line must come from these Cu 3+ sites. Given the resistiv- 
ity and specific heat evidence of a fully insulating ground 
state, they should represent 33 % of the total number of 
sites. 

It should be noted that it was not possible to single 
out the spin-singlet Cu 2+ sites which are expected to 
be present in large clusters where the local environment 
remains identical to the pristine kagome lattice as well as 
in small clusters with an even number of spins (such as 
shown in Fig. ITU]) . 

The insulating ground state of LaCu02.66 below 
^150 K is in agreement with all measurements and es- 
pecially with the high values of the resistivity and the 
absence of a linear- in-T contribution to the specific heat. 
Above 180 ±20 K, however, the resistivity, the bulk mag- 
netization and the spin-lattice relaxation of the Cul site, 
all show an activated behavior. This is reminiscent of 
semiconductors for which the thermal activation arises 
from the thermal excitations of carriers in the valence 
band. Around 180 K, there is also a crossover in the 
position and the width of the absorbtion line in electron 
spin resonance. 



B. From numerical simulations 

Wc have seen that, for a dilution level as large as 33 %, 
the identity of the original kagome lattice is seriously de- 
graded and, despite remaining patches of local kagome 
character, the system becomes mostly an inhomogencous 
pattern of magnetic clusters. Therefore, the peculiar 
properties found for weak dilution in S = 1/2 kagome 
systems (each impurity binds strong spin singlets in its 
vicinity^) are unlikely to be relevant for LaCu02.66- in- 
stead, the properties of each spin next to a Cu 3+ will pri- 
marily be determined by the morphology of the part of 
the cluster it belongs to. This does not necessarily mean 
that all memory of the kagome physics is lost. For exam- 
ple, in the S = 3/2 kagome-bilayer compound SCGO, re- 
markable properties associated with the spin liquid state 
(such as the original magnetic excitations responsible for 
a T 2 contribution to the specific heat or the peculiar dy- 
namics probed by muon spectroscop y 27 ' 37 ) are still ob- 
served at a dilution level as high as 60 %, i.e. above the 
percolation threshold p c =0.5— I 



pothesis that magnetic order is nucleated at those sites 
for which the local environment is mostly non-frustrated 
(chain-like for instance) readily justifies the partial char- 
acter of the frozen phase. This scenario of locally released 
frustration also provides a natural boost of the onset tem- 
perature of magnetic freezing which is needed to explain 
values as high as 40-50 K in LaCu02.66, that is equivalent 
to at least J/10 (see discussion on J in section III.B). 

In any case, it is likely that magnetic order in these 
frustrated compounds arises from the conjunction of dif- 
ferent factors and not from a single one. For exam- 
ple, the probable inter-triangle exchange interactions in 
YCUO2.66 represent a deviation from the simple nearest- 
neighbor kagome model which lowers the degree of 
frustration and thus favors spin ordering. This pre- 
sumably explains why the magnetic volume fraction in 
YCuC>2.66 is as high as 80 %^ much higher than 20 % in 
LaCuC>2.66j for the same dilution level. However, orbital- 
ordering and/or J arm- Teller effects might have to be con- 
sidered as well. 



VII. DISCUSSION 

A. Origin of the partial magnetic freezing 

First, it must be realized that the partial magnetic 
freezing does not arise from parasitic phases or from 
strong inhomogeneity of the sample: the former is ruled 
out by (X-ray and neutron) diffraction data- and the 
latter is inconsistent with the observation of well-defined 
NQR lines and of a well-defined structure (recall that 
oxygen dopants are not randomly located, they are well- 
ordered in Cu planes where they build the kagome lattice 
itself). Also, the well-defined oscillation seen in /iSR for 
both LaCuC>2.66 and YCuC>2.66 (with an even larger mag- 
netic fraction for the latter) speak in favor of an intrin- 
sic magnetism from macroscopic-size clusters— £ Thus, 
the frozen magnetic fraction is an intrinsic property of 

LaCu0 2 .66- 

Even if the two-dimensional, frustrated and quantum 
nature of the kagome lattice of spins 1/2 is unfavorable to 
long range ordering^, the occurrence of magnetic freezing 
in LaCuC>2.66 is not really a big surprise: inhomogeneous 
magnetism with partial fractions of frozen spins is al- 
ready found in relatively pure compounds approaching 
the perfect kagome geometry. 39 ' 40 Part of the explana- 
tion probably lies in the presence of the Dzyaloshinskii- 
Moriya interaction, which is known to drive the kagome 
antifcrromagnet towards magnetic ordering^ even for 
the quantum spins S = 1/2— - 

Dzyaloshinskii-Moriya interaction is also expected for 
the non centrosymmetric Cu-O-Cu bonds in LaCuC>2.66 
but since its magnitude is not known, it is difficult to 
assess its role in the spin freezing here. In any event, we 
expect that the destructuralization of the kagome lattice 
at 1/3-dilution should still play a dominant role. The hy- 



B. Localization of doped holes 

The properties of LaCuC>2.66 appear different from 
those of LaCu02.64- LaCu02.66 seems to have an in- 
sulating ground state with magnetic order on part of the 
sites, while LaCuC>2.64 apparently has features (such as a 
T-linear term in C p and a single NQR site with short re- 
laxation time T\) which suggest a metallic ground state 
without any magnetic order. 2 ' 14 ' 15 Because the perfect 
kagome lattice is realized only for <5 = 0.66, the small 
difference in the nominal planar-0 content of the two 
compounds implies a difference in their electronic struc- 
ture which can be much more effective than the difference 
in hole-doping in producing different ground states. Fur- 
thermore, given the lack of a well-defined planar struc- 
ture for S = 0.64, there is a much larger uncertainty on 
the oxygen content for this compound. Therefore, the 
two compounds may be significantly more different than 
what their similar formulas suggest. Nevertheless, the 
striking difference in physical behaviors is unlikely to be 
simply explained in terms of structural differences. In- 
deed, the larger disorder in LaCuC>2.64 should make the 
sample more insulating, while it is observed to be more 
metallic than LaCuC>2.66- 

A possible solution to this paradox is that the com- 
mensurate value of the oxygen concentration S = 2/3 
locks the otherwise metallic system into an insulator. A 
scenario of charge ordering at S = 2/3 is appealing for 
rationalizing a commcnsurability lock-in and for partially 
relieving frustration— 3 - It might also be compatible with 
the crossover detected near 150-180 K in the NQR and 
ESR measurements (FigJT^i). However, we only mention 
this scenario here as a possibility: in the absence of any 
direct evidence, charge order remains a speculation here 
and crossovers could as well be produced by hole local- 
ization at random locations. 
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VIII. CONCLUSION 



worth deeper theoretical and experimental investigations. 



We have found that LaCuC>2.66 has an insulating 
ground state which makes it a realization of the 1/3 
depleted kagomc lattice. We have found inhomogenous 
magnetism which numerical simulations suggest to arise 
from the heterogeneous distribution of cluster sizes and 
morphologies near the site-percolation threshold. Of 
course, a lot remains to be done to understand these 
complex delafossite compounds. Although apparently 
not displaying the metallic properties of sodium cobal- 
tates NazCoC>2 (a doped triangular lattice), this class of 
copper oxides appears as an uncharted playground for 
studying the doping of a frustrated lattice and would be 
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